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Abstract 

We show for the first time that collectively jammed disordered packings of three-dimensional 
monodisperse frictionless hard spheres can be produced and tuned using a novel numerical proto- 
col with packing density <j> as low as 0.6. This is well below the value of 0.64 associated with the 
maximally random jammed state and entirely unrelated to the ill-defined "random loose packing" 
state density. Specifically, collectively jammed packings are generated with a very narrow distri- 
bution centered at any density 4> over a wide density range <j> € [0.6, 0.74048...] with variable 
disorder. Our results support the view that there is no universal jamming point that is distin- 
guishable based on the packing density and frequency of occurence. Our jammed packings are 
mapped onto a density-order-metric plane, which provides a broader characterization of packings 
than density alone. Other packing characteristics, such as the pair correlation function, average 
contact number and fraction of rattlers are quantified and discussed. 

PACS numbers: 05.20.Jj, 45.70.-n, 61.50.Ah 
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I. INTRODUCTION 



The fundamental study of disordered jammed packings has a long history dating back at 
least to the pioneering work of Bernali^ on the so-called random-close-packing (RCP) state. 
That was traditionally thought to correspond to the densest "random" packing with a unique 
packing fraction (or density) <p ~ 0.64 for frictionless monodisperse three-dimensional (3D) 
spheres. That work spawned substantial research on disordered jammed packings.- - -^— 
About a decade ago,- it was argued that the RCP state is ill-defined for various reasons, 
including the fact that "randomness" was never quantified and the idea that a "random" 
packing can ever achieve a maximal density is not meaningful because one can infinites- 
imally increase the density of the putative RCP state, which indeed is dependent on the 
packing protocol and container boundaries, with imperceptible change in the amorphous 
pair correlation function. It was suggested that the RCP state should be replaced by the 
maximally random jammed (MRJ) state,- which is the one that minimizes a scalar order 
metric ip subject to the condition of the degree of jamming.—^ Studies of different order 
metrics for 3D frictionless spheres have consistently led to a minimum at approximately the 
same density ~ 0.64,-~- for collectively and strictly jammed packings in the plane 
(i.e., the "order map").— This consistency among the different order metrics speaks to the 
utility of the order-metric concept, even if a perfect order metric has not yet been identified. 
The fact that this jammed state is epitomized by maximal disorder has been advocated by 
other groups.— Bernal originally studied such jammed packings to understand the structure 
of liquids, but it is now known that 3D MRJ monodisperse sphere packings possess quasi- 
long-range pair correlations.— This property is markedly different from typical liquids with 
short-range interactions, which possess pair correlations decaying exponentially fast. Thus, 
MRJ sphere packings can be regarded to be prototypical glasses in that they possess the 
highest degree of disorder among all jammed packings with diverging elastic moduli.— 

The definition of the MRJ state implies that the density <f> alone is not sufficient to 
characterize jammed packings, since any packing configuration, jammed or not, is a point 
in the <j)-ip plane. This two-parameter description is but a very small subset of the rele- 
vant parameters that are necessary to fully characterize a configuration, but it nonetheless 
enables one to draw important conclusions.— The frequency of occurrence of a particular 
configuration is irrelevant insofar as the order map is concerned, i.e., the order map em- 
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phasizes a "geometric-structure" approach to analyze packings by characterizing individual 
configurations, regardless of their occurrence probability- 
One might argue that the maximum of an appropriate "entropic" metric (based on the 
frequency of occurrence of the packings) would be an ideal way to characterize the ran- 
domness of a packing and therefore the MRJ state. However, as pointed out by Ref. [j, 
a substantial hurdle to overcome in implementing such an order metric is the necessity to 
generate all possible jammed states in an unbiased fashion using a "universal" protocol in 
the large-system limit, which is an intractable problem. Even if such a universal protocol 
could be developed, the issue of what weights to assign the resulting configurations remains. 
Moreover, there are other fundamental problems with entropic measures. It is well known 
that the lack of "frustration"— in two-dimensional analogs of three-dimensional computa- 
tional and experimental protocols that lead to putative MRJ states result in packings of 
monodisperse circular disks that are highly crystalline, forming rather large triangular co- 
ordination domains. Because such highly ordered packings are the most probable outcomes 
for these typical protocols, "entropic measures" of disorder would identify these as the most 
disordered, which is clearly a misleading conclusion. 

On the other hand, an appropriate geometric-structure order metric is capable of iden- 
tifying a particular configuration not an ensemble of configurations of considerably lower 
density. For example, a jammed vacancy-diluted triangular lattice packing (and its multido- 
main variant) or a jammed packing containing many small crystalline regions and "grain" 
boundaries that is consistent with our intuitive notions of maximal disorder possesses small 
scalar order metrics. However, typical packing protocols would almost never generate such 
disordered disk configurations because of their inherent implicit bias toward undiluted crys- 



tallization. The readers are referred to Ref. 
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for further discussion. Thus we seek to 



devise order metrics that can be applied to single jammed configurations, as prescribed by 
the geometric-structure point of view. The geometric-structure approach incorporates not 
only maximally dense packings (e.g., Kepler's conjecture) and random "Bernal" packings, 
but an infinite class of other significant jammed states not previously recognized, including 
"tunneled" crystals that are putatively at the jamming threshold with <fi « 0.49365 . . .. 

Furthermore, the geometric-structure approach naturally incorporates the algorithmic 
variability of different packing protocols that leads to a diversity of density and disorder 
in jammed sphere packings.— In a typical numerical packing protocol, either the particle 
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growth or system compression leads to an increase of ^iLiiJ&iZ which makes the particle-pair 
nonoverlap constraints consume larger and larger portions of the configuration space.— Fur- 
ther increase of causes the available configuration space to fracture, generating isolated 
"islands" that each eventually collapse into final jammed states with distinct density and 
structure. Presumably, any protocol would sample the disconnected regions of the configura- 
tion space with fixed probability, leading to well-defined average of any packing characteristic 
of interest.— Unless chosen to be highly restrictive, a typical packing protocol applied to a 
system with N particles could produce a large number of geometrically distinguishable pack- 
ings with some dispersion in their and ip values. A narrowing of the distributions of the 
packing characteristics with increasing N can be expected due to operation of a central limit 
theorem. The particular values to which the distributions individually converge are protocol- 
specific, meaning that these values can be controlled by choosing appropriate protocols or 
tuning the parameters of a protocol. Indeed, jammed packings with a diversity of density 
and disorder have been produced via a variety of protocols,-^-^ 1 ^ 1 ^ including fh 0.64 
which is empirically the outcome of a considerable large number of laboratory experiments 
and numerical simulations for identical frictionless spheres. Since the jammed packings with 
« 0.64 are apparently the "most probable states", significance has been attached to the 
so-called unique J (jamming) point, i.e., ~ 0.64, which is suggested to correspond to the 
onset of collective jamming in soft-sphere systems.— However, the wide spectrum of density 
values [e.g., G (0.63, 0.74048 . . .)] that has been achieved clearly suggests that conclusions 
drawn from any particular protocol are highly specific rather than general and thus claims 
of uniqueness of packing states based on their frequency of occurrence overlook the wide 
variability of packing algorithms and the distribution of configurations that they generate. 
We will elaborate on this issue in the Conclusions and Discussion (Sec. IV). 

In this paper, we show explicitly that exploring algorithmic variability of packing pro- 
tocols can lead to a diversity of density and disorder of jammed sphere packings, which is 
consistent with the geometric-structure approach. In particular, we demonstrate that collec- 
tively jammed packings can be generated with a narrow distribution centered at any density 
over a wide range G [0.6, 0.74048. . .]. A novel sequential linear programming (SLP) 
packing algorithm^ is used to produce jammed disordered packings with as low as 0.6 
for the first time, i.e., the onset of disordered jamming occurs well below the MRJ density 
of about 0.64. The Lubachevsky-Stillinger (LS) packing algorithm^ is employed to produce 
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jammed packings with spanning continuously from that of MRJ state (0 ~ 0.64) all the 
way up to the face-centered-cubic (fee) close-packed density (0 = 0.74048 . . .).- ~— However, 
the standard LS algorithm as well as all previously used numerical protocols do not pro- 
duce disordered collectively jammed states with well below 0.64. We control the jamming 
density by tuning certain parameters of the packing protocols. The jammed packings with 
a diversity of disorder are mapped onto a density-order-metric plane. Our results strongly 
support the view that there is no universal jamming point that is distinguishable based on 
the packing density and its frequency of occurence. Other packing characteristics, such as 
the pair correlation function g 2 , average contact number Z, and fraction of rattlers f r are 
quantified and discussed. 

II. PACKING PROTOCOLS AND CONTROL PARAMETERS 

A. Lubachevsky-Stillinger Algorithm 

The LS algorithm is an event-driven molecular dynamics in which particles can grow in 
size at a certain expansion rate 7 = \dD jdt (D is the diameter of the sphere) in addition to 
their thermal motion.- Note that 7 is a key parameter that controls the jamming density.—^ 
Sufficiently small 7 enables the system almost to be in equilibrium during the densification, 
which will finally crystallize in three dimensions into a fee packing. For large 7, the system 
will quickly fall out of equilibrium and reaches a jammed state with an amorphous structure. 
Intermediate 7 will result in various degrees of partial crystallization in the system, which 
leads to a continuous range of jamming densities. It is noteworthy that a very small value 
of 7 should be used toward the jamming limit such that a true particle contact network can 
be formed for both crystalline and disordered packings.—^ We use a modified version of the 
LS alogorithm^ to generate jammed packings for > 0.64. 

B. Sequential-Linear-Programming Algorithm 

A recently devised SLP algorithm^ is used here to produce disordered jammed sphere 
packings with < 0.64. This algorithm solves an optimization problem called the adaptive 
shrinking cell (ASC) scheme:— 1 ^ jammed particle packings are generated by maximizing 
the packing density subject to interparticle nonoverlapping constraints. The optimization 
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variables include particle positions as well as shape and size of the simulation box. Starting 
from an initial configuration, a new configuration is obtained by (locally) maximizing the 
density via both individual particle motions and collective motions induced by the defor- 
mation/shrinkage of the simulation box. For spheres, the objective function and constraints 
can be linearized for a given packing configuration and the SLP method is used to solve 
the optimization problem. The final packings produced by the SLP are at least collectively 
jammed due to its incorporation of inherent collective particle motions. Details of this al- 
gorithm and its applications for generating both disordered and maximally dense packings 
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in high dimensional Euclidean space are given in Ref. 

By removing spheres from the fee packing and its stacking variants without destroying the 
rigidity of the contact network, Torquato and Stillinger— have constructed strictly jammed 
"tunneled crystal" packings with <fi as low as 0.49365 .... Similar removal procedures can be 
applied to MRJ packings. However the contact network of the remaining packing is generally 
not rigid anymore. The removed particles are required to have a large coordination number 
and to be mutually separated by at least a few sphere diameters. Compressing the remaining 
packing using the SLP algorithm leads to a re-jammed configuration, with a minimum 
degree of structural relaxation (i.e., significant particle reconfigurations are localized around 
the cavities).— The re-jammed configurations generally possess < 0.64 and the removal- 
compression procedure can be repeated several times until a lower limit of <j) is reached. 
The control parameters of the SLP algorithm include the initial MRJ configurations and the 
number of removed particles. 



III. JAMMED PACKINGS WITH VARIABLE DENSITY AND DISORDER 

A. Histograms for Jammed Packings with (f> > 0.64 

We employ the LS algorithm to generate a large number of jammed packings of monodis- 
perse spheres in three dimensions for <fi > 0.64 with N = 250, 500, 1000, 2500, 5000, and 
10000, and a wide range of initial expansion rate 7 € [10~ 2 , 10~ 6 ]. Results are verified to be 
at least collectively jammed,— using either a linear programming protocol or by monitoring 
the instantaneous pressure of the systems for long time periods;— and rattlers are included 
to compute the reported densities. By tuning 7, the density at which the systems jam can 
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be controlled. Figure Hfb), (c) and (d) contrast distributions of for two distinctly different 
system sizes (N = 250 and 2500) converging onto ~ 0.64, 0.68 and 0.72, respectively. (The 
case of ~ 0.60 shown in Fig. 1(a) will be discussed separately below.) It is clear that as N 
increases the distributions narrow for all three mean density values. We also find that such 
narrowing becomes even more significant for N = 5000 and 10000. Specifically, Fig. [2] shows 
how the standard deviation a of the density distribution varies with system size N for the 
case in which the mean « 0.66. Observe that a is a monotonically decreasing function of 
N and a approximately scales N^ 1 ^ 2 for large N. A similar narrowing of the distribution 
occurs for « 0.60 shown in Fig. 1(a), and this case will be discussed separately below. 
Thus, one can expect that in the "thermodynamic" limit (i.e., N — > oo) the jamming density 
will converge to a well-defined value anywhere over the interval £ [0.60, 0.74048 . . .]. These 
results imply that any inclination to select a specific value as uniquely significant (e.g., 
pa 0.64) is primarily based on inadequate sampling of the full range of algorithmic richness 
and diversity that is available at least in the underlying mathematical theory of jamming. 

B. Histograms for Jammed Packings with < 0.64 

We produce the majority of jammed packings with < 0.64 using the SLP algorithm. 
MRJ packings generated via the LS algorithm are used as initial configurations. Each time 
approximately f s = 0.1% — 2.5% of the spheres are removed from the initial packing and 
the remaining spheres are compressed to a jammed state using the SLP algorithm. Such a 
procedure is repeated n r = 5 — 10 times before a lower limit on is reached. The fraction 
f s of removed spheres is decreased as the limit is approached. We stress that our packings 
with < 0.64 are not so-called random loose packings,— which are not even collectively 
jammed.— A few packings with pa 0.62 are generated using the LS algorithm with open 
simple-cubic lattice packings as initial configurations.- 

Figure Ufa) contrasts distributions for two distinctly different system sizes (N = 216 
and 2235) converging onto ~ 0.60. It can be seen that as N increases the distributions 
narrow. Such narrowing is also observed for other convergence densities within the interval 
[0.6, 0.64]. We note that it is very difficult to produce jammed packings with significantly 
lower than 0.6 using the SLP algorithm (though removing the rattlers in the packing results 
in a slightly lower density ~ 0.595), while it has been rigorously shown that strictly jammed 
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sphere packings (e.g., tunneled crystals and the associated stacking variants) can possess 



and MRJ packings can be constructed. In particular, stackings of layers of honeycomb-lattice 
packings of spheres stabilized with triangular-lattice layers on top and bottom are inserted 
into the MRJ packings. These "layered" packings are then compressed to jamming using 
the SLP algorithm. This construction enables one to obtain jammed packings with variable 
disorder within the density range 6 (0.49,0.64). 

C. Order Metrics and Other Packing Characteristics for Jammed Packing with 

€ [0.49, 0.74] 

It has been established that a variety of different useful order metrics are positively 
correlated^ - — and hence any one of them can be used to characterize the packings. To 
quantify the order of the packings, we compute the translational order metric T,- defined as 



where rij is the average occupation number for the shell i centered at a distance from a 
reference sphere that equals the zth nearest-neighbor separation for the open FCC lattice at 
that density and N s is the total number of shells for the summation (N s = 45 is used here 
for N ~ 2500); rif eal and nf cc are the corresponding shell occupation numbers for an ideal 
gas (spatially uncorrelated spheres) and the open FCC lattice. For a completely disordered 
system (e.g., a Poisson distribution of points) T = 0, whereas 1=1 for the FCC lattice. 

Figure [3] shows the 0-1 plane on which representative jammed packings with N ~ 2500 
and G [0.49, 0.74] are mapped. The packings generated using the LS algorithm are 
shown as black circles. Note that for these packings each is associated with a range of X 
values. Moreover, increasing from the MRJ state can be achieved at the cost of decreasing 
the degree of disorder, as pointed out in Ref. Q. Further increase of is associated with 
partial crystallization of the packings, as indicated by the sharp peaks of the pair correlation 
function gi of the packings (see Fig. Hj). The fraction of "rattlers" (i.e., locally unjammed 
individual particles that can move freely within cages of jammed neighbors) decreases (from 
~ 2.8% with pa 0.64 to 0% with = 0.74048. . .) and the average contact number per 
particle Z increases (from 6 to 12) as the density increases due to the (partial) crystallization 



as low as :y -| 2E = 0.49365 . . ..— However, jammed packings that combine the tunneled crystals 
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of the packings, and we have Z w 6.00, 6.45 and 7.96 respectively for densities ~ 0.64, 0.68 
and 0.72 (rattlers are excluded when computing Z). 

Several representative packings obtained via the SLP algorithm are also mapped onto the 
(f>-% plane (red squares in Fig. [3]).— It can be seen that jamming at <ft < 0.64 is necessarily 
associated with an increase in the degree of order, which is also indicated by the increase 
of average contact number per particle (i.e., Z m 6.37 for <fi « 0.6) and the decrease of 
the fraction of rattlers (i.e., ~ 1.1% for <p « 0.6). The pair correlation function g2 of a 
representative packing configuration (with N = 2235) is shown in Fig. HJ 

The "tunneled crystal" packings are mapped onto the <j) — X plane (green up-triangles) . 
The FCC tunneled crystal possesses the highest translational order metric X, the "disor- 
dered" Barlow tunneled crystal (a random stacking of layers of honeycomb-lattice sphere 
packings) possesses the lowest X, and the "zig-zag" tunneled crystal (analog of the hexagonal- 
close sphere packing)^ is in between. The dashed blue lines show the spectrum of packings 
generated by randomly filling the vacancies in the corresponding "tunneled crystal" packings, 
leading to perfect FCC, HCP and "disordered" Barlow packings at the maximal density- 
Several representative layered packings are also mapped onto the 0-X plane (purple dia- 
monds) in Fig. [31 The simple-cubic (<fr = 0.5235 . . ., blue left-triangle) and body-centered- 
cubic (</> = 0.6801 . . ., orange right-triangle) packings are also shown in Fig. [3j Note the SC 
and BCC packings are not collectively jammed. 



IV. CONCLUSIONS AND DISCUSSION 



In summary, we have shown that jammed sphere packings can be produced with a narrow 
distribution centered at any <p over a wide range <fi 6 [0.6, 0.74048...] with a diversity 
of disorder, by exploring the algorithmic variability of packing protocols. This suggests 
that any temptation to select a specific <fi value as uniquely significant or universal (e.g., 
(f) ps 0.64) based on its frequence of occurrence is primarily due to an inadequate sampling 
of the full range of algorithmic richness and variability of packing protocols. Our packings 
are characterized by a "geometric-structure" approach, i.e., they are mapped onto the 0-X 
plane; and we have shown that moving away from the MRJ density in both directions (i.e., 
lower or higher <p) leads to a higher degree of order and a larger average contact number. 

It is claimed that the protocol (e.g., the athermal relaxation method used in Ref. llOl ) 
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in which the jammed states of soft spheres are weighted by the volume of their basins of 
attraction has a clear interpretation in terms of the energy landscape, e.g., an ensemble of 
jammed states is considered. In fact, the conceptual framework of such a protocol (e.g., 
the energy landscape picture and the ensemble of inherent structures) was introduced by 
one of us over 40 years ago.— 7 — Although this approach is well-defined and theoretically 
possible, it is inevitable that all protocols sample the energy landscape in their own biased 
fashion. Therefore, there is no compelling mathematical criterion for selecting one protocol 
over the others. For example, from completely random initial configurations (i.e., Poisson 
distribution of points), it is found the jamming density sharply peaks at <j) ~ 0.64.— However, 
it was shown nearly a decade ago that a much wider density range can be achieved for both 
monodisperse and polydisperse sphere packings using the Lubachevsky-Stillinger packing 
algorithm;- - -^ and recently similar results for polydisperse spheres have been obtained^ 
using the same athermal protocol as in Ref. |l0|. More recently, we have devised a novel 



(athermal) linear-programming packing protocol that enables one to obtain a spectrum 
of inherent structures ranging from disordered jammed packings up to the maximal density 
packings starting from random initial configurations.— All of these results suggest that there 
is no "universal" protocol that can generate all possible jammed states in an unbiased 
fashion. 

Moreover, it is not clear how the jammed states should be weighted in an ensemble. It has 
been suggested that the volume of the basins of attraction for systems starting from Poisson 
distributions of interacting points^ should be used for the weighting. However, in reality 
virtually no jammed systems are experimentally produced using a Poisson distribution as 
an initial configuration. In addition, the jammed packings produced either numerically or 
experimentally always contain a small but different fraction of rattlers. Strictly speaking, 
these jammed states are not single points on the energy landscape but bounded regions with 
dimensions equal to the number of degrees of freedom of the rattlers. The dimension of 
the jamming basins are different, and therefore it is ambiguous to consider their volumes 
for weighting and to characterize them on the same footing. Ideally, each jamming basin 
should be characterized individually as emphasized in the "geometric-structure" approach. 
Removing the rattlers will not affect the jamming nature of the packings, but leads to the 
lower jamming density. This further adds ambiguity to a density-alone characterization of 
jammed packings, which is employed in the "ensemble" approach. 
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We also would like to note that as exercised the "ensemble" approach has focused on 
protocols that disproportionately generate disordered packings, presumably because such 
packings have a high frequency of occurrence in typical experiments or simulations. The 
specific protocols employed discriminate against the ordered structures such as the FCC 
packing (and its stacking variants) and the tunneled crystals^ and result in an inappropri- 
ate narrow vision of the set of all possible collectively jammed packings. Thus, the physical 
relevance of jammed packings should not be determined based on their frequency of occur- 
rence in experiments or simulations, which again are protocol dependent. 

From our point of view, the "ensemble" and "geometric-structure" approaches do not con- 
flict with each other but rather are complementary. For example, the "geometric-structure" 
approach characterizes individual packing configurations drawn from well-defined ensembles. 
However, in our experience, the "geometric-structure" approach can always provide nontriv- 
ial solutions when the "ensemble" approach breaks down, such as the possible identification 
of disordered jammed two-dimensional packings discussed in the Introduction. Therefore, 
we believe that the "geometric-structure" approach, when utilized in the broad context of a 
full set of available protocols, has a capability to discover, integrate and characterize packing 
structures that go well beyond the disordered set. We emphasize that all of the aforemen- 
tioned issues have been addressed in the literature.-~-^22, In this paper, we have examined 
the applications of the "geometric-structure" approach in a broader context by amplifying 
and quantifying some of these issues. 

An interesting open question that naturally arises is whether <fi m 0.60 is the lower 
limit on the density of jammed amorphous sphere packings? We note that the distribution 
of density for N ~ 2500 shown in Fig. Ufa) can be fitted with a Gaussian with mean 
4> = 0.602 and standard deviation a = 0.002. This strongly indicates that there is a 
small but finite probability of finding jammed packings with even lower densities. Given 
enough number of trials, such packings could be obtained in principle. It is noteworthy 
that the density <p = -^p = 0.49365 . . . associated with the "tunneled crystals" is likely to 
be the threshold (i.e., lowest possible) density for strictly jammed sphere packings in three 
dimensions.— However, neither the LS nor the SLP algorithms as normally implemented 
is able to produce such packings, presumably because both algorithms tend to densify the 
packing and jamming is a consequence of their "compression" nature. Lower <fi is attainable 
via the SLP algorithm because jamming is achieved through local reconfigurations. However, 
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it is currently not designed to find highly unsaturated jammed packings, such as the tunneled 
crystals. Therefore, we see no reason that jammed amorphous sphere packings with even 
lower density cannot be produced via carefully designed protocols. 

A limitation of current protocols designed to produce jammed packings is that they in- 
evitably lead to packings with high densities. It is highly desirable to devise protocols that 
explicitly take into account the requirement of jamming as well as other packing charac- 
teristics (e.g., <p and Z). One possible approach is to delineate the conditions for jamming 
(and other characteristics) in a quantitative way and to include them as constraints of an 
optimization problem. We have suggested possible solutions to the problem of producing 
low-density jammed amorphous packings in Ref. 
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In future research, we will focus on the 



development of such packing protocols, which is a highly nontrivial challenge. 
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FIG. 1: Histograms of jammed sphere packings that are centered around different mean densities. 
Panel (a) shows packings generated using the SLP algorithm with eft ~ 0.6. Panels (b), (c) and (d) 
show packings generated using the LS algorithm with <f> ~ 0.64, 0.68 and 0.72, respectively. The 
distributions become narrower as the system size increases. 

FIG. 2: The standard deviation a of the density distribution with mean (f> ~ 0.66 as a function of 
system size N . System sizes N = 500, 1000, 2500, 5000 and 10000 are used here. The linear fit for 
lner vs. In TV" gives a slope k ~ —0.533, which indicates that a approximately scales as N~ 1 ' 2 for 
large N. Such a a-N relation is also observed for density distributions with the other mean values 
studies here. 

FIG. 3: Order map for jammed sphere packings with N ~ 2500: translational order metric T versus 
packing density (p. The dashed blue lines (which are consistent with the qualitative trends indicated 



in the order maps discussed in Ref. Il3l ) show the spectrum of packings generated by randomly 
filling the vacancies in the corresponding tunneled crystal pekings, which leads to perfect FCC, 
HCP and Barlow packings at the maximal density. The simple-cubic (SC) and body-centered-cubic 
(BCC) packings (not jammed) are also shown. 

FIG. 4: Pair correlation function 52 and average contact number Z of representative jammed sphere 
packings at different densities, where D is the sphere diameter. 
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